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A SYSTEM FOR COMPACTION AMD RECONSTRUCTION OF WAVELET DATA 

FIELD OF THE INVENTION 

The present inveniion lebtes to data processing systems generally and to 
systems utilizing wavelet transforms in partictilar. 

BACKGROUND OF THE INVENTION 

Data proccssmg typically involves transforming the data 6om one space lo 
another. The best known transformation is the Fomier transformation which transforms a 
signal from a time space to a frequency space, indicating thereby the frequency components 
of the input signal. 

Recently, wavelet transforms have been developed where each wavelet 
transform has a different curving "wavelet" shape and the result of the transformation 
indicates to what extent the input signal, a digital signal, is formed of a combination of the 
different wavelets. The article "Wavelets and Signal Processing", by Olivier Rioul and Martin 
Vetlerli, published in TFEF SP" Magazine. October 1991, pp. 14 - 38, describes prior art 
wavelet transform operations and is incorporated herein by reference. 

Typically, wavelet transforms are implemented in a multi-resolution marmer. 
The prior art multi-resolution method and exemplary output, for a one-dimensional signal, are 
illustrated in Figs. 1 A and 2, respectively, to which reference is now made. 

Shown in Fig. lA are three "wavelet transform units" lOa - 10c. Each wavelet 
uansform unit lOa - 10c is formed of its respective wavelet low pass filler (WLPF) 12a - 12c, 
wavelet high pass filter (WHPF) 14a - 14c and samplers 16a - 16c. The wavelet low and 
high pass filters 12a - 12c and 14a - 14c break the signal into its high and low wavelet 
frequency components and the samplers 16a - 16c sample the output of each filter, thereby 
changing the resolution of the signal. For one-dimensional input signals, the samplers 16a - 
16c typically take every second datapoint, thereby reducing the number of points by two with 
each level. Other sampling rates can also be implemented. 

Each transform unit lOa - lOc operates on a different resolution signal. 
TTierefore, the elements of each transform unit are labeled to indicate the resolution level or 
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iteration. 1 - 3, at which they opeiate. wheiein the highest level has the lowest lesolution 
Oeast number of points). The fiist transform trnit 10a operates on the input signal and the 
other transform units 10b - lOc operate on the low ftequency output signals of the transfonn 
unit 10a - 10b, respectively, of the level previous to it. The signals of Fig. 1 are labeled to 
indicate the level and whether they are the low (L) or high (H) frequency output of the level. 

The output of a wavelet uansform operation comprises the foUowing signals, 
the high frequency output signals HI - H3 and the low frequency output signal L3 of the last 
level. The low frequency output signal of the highest level, in this example signal L3, is 

herein called the "core" signal. 

An exemplary input sipal, labeled 18. and its conesponding three high 
frequency output signals are shown in Fig. 2. As can be seen, the resolution is reduced at 
each level. Furthermore, the activity of the lower resolution signals is near at the inflection 

point of the input signal 18. 

If the input signal 18 is an image or some other two-dimensional signal, the 
15 wavelet low and high pass filters can either be separable or non-separable. The operations 
for separable and non-separable filieis are shown in Figs. IB and IC, respectively, to which 

reference is now made. 

The wavelet transform units 10 for separable filtering include low and high 
pass wavelet fillets 11a and lib (Fig. IB) and samplers 13a for rows and low and high pass 

20 wavelet fillers 1 Ic and 1 Id and samplers 13b for columns. The rows of the input signal fust 
pass through the low and high pass wavelet filters 11a and lib and are sampled by the 
sampleis 13a, thereby reducing the number of pixels in each row by two. The columns of 
the row fjiered signals then pass through the low and high pass wavelet filters 1 Ic and lid 
and are sampled by the samplers 13b, thereby reducing the number of pixels in the columns 

25 by two. Alternatively, the columns can be operated on fust, followed by the rows. 

The separable filtering method produces three high frequency signals Ha, Hb, 
He and one low frequency signal L at each level. Thus, the fust level produces the signals 
Hla, Hlb, Hlc and LI. Each signal has one-fourth the number of pixels as the input signal. 
The low frequency signal LI is then provided as the input sipal to the next wavelet 

30 transform unit 10. 

The wavelet transform units 10 for non-separable filtering operate differently 

for even- and odd-numbered levels, as illusirated in Fig. IC. For odd levels, the input signal 
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is operated on by non-separable low and high pass wavelet filters 15a and 15b, respectively, 
and quincunx odd samplers 17a. Odd samplers 17a perform "quincunx" subsampling which 
removes every even daiapoint in odd numbered rows and every odd datapoint in even 
numbered rows, leaving a checkerboard-like output image, illustrated in part in Fig. lOB. In 

5 order to maintain the image, the datapoints are not actually removed; instead, they just arc 
not operated upon by the foUowing even level. 

For even levels, the wavelet transform unit 10 operates on the checkerboard- 
like image with non-separable low and high pass wavelet filters 15c and 15d, respectively, 
and quincunx even samplers 17b. Even samplers 17b also perform quincunx sampling on the 

10 checkerboard-like image, resulting in a more normal looking image, shown in part in Fig. 
lOA, which is the input to the odd-level wavelet transform unit The image produced by even 
samplers 17b has one-half the number of elements as the input image to the odd samplers 
17a. 

The article "Predictive Interscale Image Coding Using Vector Quantization", 
15 by M. Antonini et al and published in 5;ipna1 Piocessing VtTheories and A ppb'calions, L. 
Torres, E. Masgrau and M.A. Lagunas, editors, Elsevier Science Publishers, B.V. 1990, pp. 
1091 - 1094, describes an image compression method using biorlhogonal wavelet transforms 
and an interscale prediction scheme. The prediction scheme predicts the location and 
amplitude of edges in an image at one resolution based on their locations and amplitudes al 

20 the next highest resolution. 

The article "Image Coding using Lattice Vector Quantization of Wavelet 
Coefficients", by M. Antonini et al., ICASSP 91. pp. 2273 - 2276 describes a compression 
method for wavelet transform data. 

U.S. Patent 5,014,134 to Lawion el al. describes an image compression and 
25 decompression system which utilizes wavelet transforms. 
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SUMMARY OF THE PRESENT INVENTION 



It is an object of the present invention to provide a system and method of 
compaction and leconstmction of wavelet transform data which preserves areas of high 

5 information content in an inpnt signal. The input signal can be a one-dimensional signal, a 
two-dimensional image or a many dimensional signal. 

There is therefore provided, in accordance with a prefened embodiment of the 
pxesent invention, a system for compaction and reconstruction of wavelet transform data of 
an input signal where the wavelet transform data is formed of a pluiaUty of resolution levels. 

10 The system includes a compaction unit and a convener. The compaction unit is operative on 
forests of trees, wherein each Uee has branches and each branch relates one datapoint in a low 
resolution level with the entirety of its descendant datapoints in higher resolution levels. The 
compaction unit prunes branches from the trees and trees from the forest to remove all but 
areas of high information content in the input signal, thereby to produce pnmed trees. The 

15 convener convens the primed trees to reconstnacted wavelet transform data. 

Additionally, in accordance with a prefened embodiment of the present 
invention, the compaction unit includes a) a non-null pnming unit for removing branches of 
the trees having datapoints thereon whose wavelet uansfonn values are below predetermined 
threshold levels, producing thereby non-nuU pnmed uees and b) a zerc^crossed pnming unit 

20 for detennining locations of transients in the input signal and for creating zero-crossed pnmed 
uees having therein only datapoints of the non-null pnined trees which are at or near the 

locations of uansients. 

Moreover, in accordance with a prefened embodiment of the present invention, 
the datapoints of the zero-crossed pnmed uees fonn nodes and the compaction unit further 
25 includes apparatus for associating datapoints of the non-Dull pnined Uees with nodes of said 
zero-crossed pnined uees, wherein each datapoint can be within a local area of many nodes 

but is associated with only one of them. 

Funher, in accordance with a prefened embodiment of the present invention, 
the system includes a nullifier, operative prior to the compaction unit, for nulling those 
30 datapoints of the wavelet uansfonn data which are below the predetennined threshold levels. 

StiD funher, in accordance with a prefened embodiment of the present 
invention, the zero-crossed pnming unit includes apparatus for detennining locations of zero- 
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ciossings in the wavelet transfonn data. 

Finally, the method perfonns the operations of the elements of the system. 
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BRIEF DESCRIPTION OF THE DRAWINGS 



The preseDl invention will be undeisiood and appreciated more fiJly from the 
foDowing detailed description taken in conjunction with the drawings in which: 
5 Pig. lA is a block diagram illustration of a prior multi-resolution wavelet 

transform operation; 

Fig. IB is a block diagram illustration of one section of a prior art wavelet 
transform operation for separable filtering; 

Fig. IC is a block diagram illustration of two sections of a prior art wavelet 
10 transfom operation for non-sepaiable filtering; 

Fig. 2 is a graphical illustration of an input signal and versions of the signal 
after wavelet uansformation having high frequency information therein; 

Fig. 3 is a block diagram illustration of a wavelet transform compaction and 
reconstruction system, constructed and operative in accordance with a preferred embodiment 

15 of the present invention; 

Fig. 4 is a graphical illusuation of the parent-child relationships between the 
signals of Fig. 2, useful in understanding the operation of the present invention; 

Figs. 5A and 5B are schematic illustrations of full and pruned trees, useful in 
understanding the operation of the system of Fig. 3; 
20 Fig. 6 is a block diagram iUusiraiion of a method of nullifying some wavelet 

transform data, useful in the system of Fig. 3; 

Fig. 7 is a flow chart illustration of the operation of a compaction unit forming 

part of the system of Fig. 3; 

Fig. 8 is a flow chart illustration of a method of producing a zero-crossed tree, 

25 forming part of the operations of the compaction unit of Fig. 7; 

Fig. 9 is a flow chart illustration of a method of computing zero-crossing 
locations forming part of the method of Fig. 8; 

Fig. lOA is a schematic illusuation of a local area wdthin a two-dimensional 

image; 

30 Fig. lOB is a schematic iDusuation of a local area for an odd level wherein 

quincunx subsampling was performed; and 

Fig. 11 is a flow chart illustration of a method of producing a fully pruned tree 
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foimiBg part of the method of Fig. 7. 



DETAILED DESCRIPTION OF A PREFERRED EMBODIMENT 



The present iDvention will be described hereinbelow in the context of one- and 
5 two-diniensional input signals, it being understood that input signals of other orders of 
magnitude are also included vvithin the scope of the present invention. 

Reference is now made to Fig. 3 which illustrates a wavelet tiansfonn 
compaction and reconstruction system, constructed and operative in accordance with a 
prefened embodiment of the present invention. The system of the present invention typically 
10 comprises a compaction system 20 for compacting wavelet transform signals of a digital input 
signal and a leconstruction system 22 for reconstructing the input signal from the compacted 

wavelet transform signals. 

The compaction system 20 typicaUy comprises a wavelet transformer 24, a 
Dullifier 26 for zeroing out values of the wavelet transform data W which are below a 
15 predetennined human-defined threshold and a compaction unit 28 which compacts the output 

W of the nuUifier 26. 

The wavelet transfomcr 24 comprises a predetermined plurality of wavelet 
transform units 10 combined as described in the Background of the Present Invention. For 
most of the examples herein, wavelet transformer 24 comprises three wavelet uansform units 
20 10, however it will be understood that other numbers of units 10, typicaUy in the range of 

three to eight, are suitable. 

The v^'avelet transform units 10 can be separable or non-separable. It is noted 

that the type of wavelet transform utilized (either separable or non-separable) affects the 

percentage of compaction to be achieved, as will be described hereinbelow. 
75 As is known in the art, the wavelet transformer 24 produces two types of 

sipals. the core signal and the detailed signal. The core signal is the low frequency signal, 

labeled L3 in Fig. 1, which has a low resolution. The detailed signal comprises the high 

frequency signals HI - H3 each having a different resolution. 

If separable filtering was utilized on a two-dimensional input image, each high 
30 frequency signal HI - H3 is formed of three signals. For the remainder of this discussion, 

the term HI will be utilized and it will refer to the output of both separable and non-separable 

filtering unless otherwise specified. 
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m compaction imit 28 l^ically detennine the areas of "high information 
content" or of detail in the nullified wavelet tiansfom data W*. These areas ean be defined 
as those areas which have transients in them and which exist in most or aU of the high 
frequency signals of the nullified wavelet data W°. 

In accordance with the present invention and as illustrated in Figs. 4, 5A and 
5B to which reference is now briefiy made, the compacUon unit 28 organizes the detailed 
signal as a "forest" of "trees" connecting the datapoints of the higher resolution signals (HI 
and H2) with their conesponding datapoints in the lowest resolution signal (H3). The 
datapoints of the lowest resolution level form the "roots" of the trees and therefore, the lowest 
resolution signal H3 is herein denoted the "root signal". There is no more than one tree for 

each daiapoint in the root signal. 

The compaction unit 28 then "prunes" the trees to remove datapoints having 
Uttle information content. The compaction unit 28 can also remove trees from the forest if 

there is no information in a tree. 

Figs. 5A and 5B illustrate two full Uees 30 and two pruned trees 32 for two 
datapoints (0 and (i+ 1) of the root signal of a one-dimensional signal. In this example, there 
are five levels and thus, five high frequency signals, HI - H5, are illustrated. In this example, 
each wavelet transform unit 10 reduced the resolution of the signal by 2. Therefore, each 
daiapoint at one resolution level conesponds to two datapoints in the next highest resolution 
20 level. 

It is noted that each tree 30 and 32 comprises a multiplicity of branches 34 
connecting each daiapoint or node 35 of a signal with its conesponding datapoints of the 
signal having the next highest resolution. It is also noted that there are no branches 34 of the 
pruned trees 32 which are not comiecied to the structure. Thus, analogous to pruning a 

25 physical tree, when one branch 34 is removed, all of the branches 34 comiected to it in the 
higher resolution signals must also be removed. 

Refening back to Fig. 3. the compaction unit 28 produces a forest of pruned 
trees 32, at most one tree for each datapoint in the root signal. The forest of pruned trees 32 
contains' generaDy all the essential information in the input signal, but in a more compact 

30 format than the output of the wavelet transform. TTie compact forest of pruned trees 32 can 
then be transmitted as desired and as illustrated in Fig. 3. For the purposes of the present 
discussion, "Uansmission" indicates transmission between any two elements, whether through 
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commmiicalion devices or by storing on a storage device. 

The leconstruction system 22 of the present invention typicaUy comprises a 
convener 40 for converting the pr^ed forest back to a wavelet transform representation W 
and an inverse ^vavekt transformer 42 for performing the inverse wavelet transform operation, 
5 thereby to reconstruct the input signal from the wavelet transform representation W, 
Because the compaction system 20 maintains generaUy the essentia] details of the input signal, 
the reconstructed signal suongly resembles the input signal. 

It vvill be appreciated that the system of the present invention can be utilized 
in any situation where the input signal must be reduced in size while the information content 
10 within it be maintained. For example, if the input signal is a high quality digital image, 
comprised of more than one-quarter million picture elements or pi^tels. it can be compacted 
with the present system for storage and reconstructed for operation. Alternatively, it can be 
transmitted to a second location after compaction and it can be reconstructed at the second 
lcM:alion- 

15 The compacted form of the wavelet transform lends itself easily to further 

compression, as foUows: Tht core signal may be further compressed by known Discrete 
Cosine Transform (DCT) methods; the fuU-pruned forest topology may be encoded by 
standard tree encoding methods; and the daiapoints of each node of the fuDy-pruned foresl 
may be considered as a vector and thus, vector quantization may be applied to compress 

20 them. 

Reference is now made to Fig. 6 which illustrates the operations of nullifier 
26. The equations performed by the nullifier 26, for a two dimensional input signal (e.g. an 
image) are listed hereinbelow in the section labeled "Nullifier Equations". The index f 
indicates the resolution level and (i,j) are spatial indices. The nullifier 26 determines a 
25 threshold level for nullifying, based on the type of input signal and on its particular features. 
Any values of the wavelet data less than the threshold value are nullified (e.g. set to 0). 

Because the wavelet transform data W is mulii-resolutional and the values in 
a signal at one resolution level are not equivalent to those in a signal at another resolution 
leve*! due to human response mechanisms to the specific type of input signal (audio, video, 
30 etc), the data W, in step 52, has to be normalized according to human responses. 

The human response function R(«,iJ) depends on the resolution level f and also 
on the position of each datapoint (y) within the image. For the type of wavelet transforms 
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described hereinabove, an exemplaiy human response function R(?,ij) is; 
R({,iJ) = R(0 = 0.95*2.2 



-«2 (1) 



Prior lo normalization, in step 50. the absolute values of the data are 
detennined, producing signal W,. The output of step 52 is a signal W,. 

After normalization, the data is quantized (step 54), with a quantization level 
Q, to reduce the number of possible values, and then clipped (step 56) to maintain the number 
of possible values to a fmite and reasonable amount, such as 1024. 

It will be appreciated that steps 50 - 56 can be performed at once, if desired. 
A histogram of the resultant normalized, quantized data from aU of the high 
frequency signals, is produced in step 58. In step 60 an overaU threshold level T„ is 
deieimined, based on the desired percentage, such as 80%. of the histogram to be nullified. 
The overall threshold level is the value of the histogram for which the desired percentage of 
15 the data has values less than or equal lo it 

The overaU threshold level cannot be utilized for the threshold operation (step 
64) since its value is a normalized one. Therefore, in step 62. the threshold is "denormalized" 
(i.e. the inverse of the normalization operation) to provide threshold values T(f,ij) for each 

resolution level f and across the image. 
20 In step 64, the threshold operation is performed wherein each datapoint in each 

high frequency signal HI - H3 of the wavelet transform data W is compared to its 
conespoDding threshold and those daiapcints which are less in absolute value than the 

conesponding threshold are nullified. 

Reference is now made to Fig. 7 which illustrates the operations of compaction 
25 unit 28. Compaction unit 28 performs two types of pruning, a first pruning (step 70) to 
remove those branches which have null data on them, and a second pruning (step 72) to 
remove those branches which do not have transient data. Since transients are typicaUy more 
than one datapoint wide, compaction unit 28 then enlarges (step 74) the nodes of the resultant 
pruned tree to include at least some of the non-nuU datapoints near the datapoint of the node. 

The pruning steps 70 - 74 can be performed by creating lists of the nodes 
which belong in each pruned tree. In accordance with an alternative embodiment of the 
present invention, each datapoint can have many attributes and pruning is peiformed by 
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updBiiBg the aitributes of each node which belongs in the pnmed tree. The steps 70 - 74 
be described as though new pnmed tiees aie produced, it being understood that either method 
can be performed. 

Prior to discussing the specific details of each operation of compaction unit 28, 
we define a parent node and its child nodes of the level below it For a one-dimensional 
signal, a parent node W^CU) has the foUowing children nodes: W<'(M,2i+l). V(M,2i). 

For a two-dimensional signal upon which separable filtering was performed, 
a parent node W<'(P,d,ij) has the following child nodes: W°(M,d,2i+l,2j+l), W0(^ 
l,d,2i-^1.2j), W<'(M,d,2i,2j>l) and W«(M,d,2i,2j). where d indicates the type of filtering 
plr'formed. For example, d = 1 indicates that the rows were low pass filtered and the 
columns were high pass filtered, d = 2 indicates that both rows and columns were high pass 
filtered and d = 3 indicates that the jo^ were high pass filtered and the columns were low 
pass filtered. 

For a two-dimensional signal upon which quincunx subsampling was 
performed, a parent node W«(f.ij) of an even-numbered level has the following child nodes 
of the odd-numbered level below it: W<'(M,2i,2j), W«(M,2i-M,2j>l). A parent node W°(fi,i.j) 
of an odd level has the foUowing child nodes of the even-numbered level below it: W°(l- 1 ,i j) 
and W°(M,i+lj) for j even and W^Cf-Lij) and W^Cf-Li-Lj) for j odd. It is noted that in the 

odd levels, if i+j is odd, the datapoint docs not exist. 
20 Similar relationships can be built for input signals of higher dimensions. 

b step 70, the compaction unit 28 reviews the entirety of nodes in the forest 
of trees. For each node which is null and whose descendants do not have significant 
information in ibem, where "significance" is defined hereinbclow in equation (2), the node 
and its descendants are removed from the tree to wWch it belongs. TypicaUy, the compaction 
25 unit 28 begins the review with the root signal H3. 

The significance S of the data in the descendants of the node is determined as 

follows: 
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30 



S - r W<'(*)R(*) 

where (*) indicates indices and denotes a descendant of the node 35, W<'(*) is the value of 
the descendant, Rr) is the value of the normalization function for the descendant, and the 
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sum 



is over all the desceDdanls of the node. S has to be above a piedetennined threshold 
value in order for the node to be kept within the pruned uee structure. 

For those nodes 35 which have nuD values, the compaction unit 28 determines 
whether or not the branches 34 attached to the node have significant data in them. If so, it 
5 saves the node as part of the pruned tree structure. 

The compaction unit 28 repeats the above determination for each of the 
datapoints in the root signal and for each of the high frequency signals. The result is a forest 

of "non-nuH" pruned trees. 

In step 72, the compaction unit 28 reviews the non-null pruned tree and the 
10 wavelet transformation W to deieimine the areas conesponding to "transients". A transient 
is herein defined as the location where the wavelet transformation value changes from positive 
to negative or vise-versa; in other words, a transient is the location of a "zero-crossing". By 
choosing the appropriate wavelet bases for the wavelet uansformation operation (performed 
by wavelet transformer 24 (Fig. 3)), the zero crossings can be forced to be at locations of 
15 maximum transitions in the input signal. 

Selection of appropriate wavelet bases is discussed in the foUowing anicle 
which is incorporated herein by reference: Stephane MaUet, "Zero-Crossings of a Wavelet 
Tr.ncfnm,", TFFF Transactions nn Informaiion Theory Vol. 37, No. 4, July 1991, pp. 1020 - 
1033. 

20 The operations performed are provided in the flow charts of Figs. 8 and 9, to 

which reference is now made. Fig. 8 iTustrates the operations performed to produce the zero- 
crossed pruned trees and Fig. 9 illustrates the operations performed to determine where the 

zero-crossings occur. 

In step 80 of Fig. 8, which is detaUed in Fig. 9, the compaction unit 28 
25 determines which nodes 35 of the non-null pruned forest are part of transients. The location 
of the transients are determined through an iteraUve process operating on the non-nullified 

wavelet uansforra W. ' 

The iteration process shown in Fig. 9 operates on a "local area" or 
neighborhood in W of each node 35 b the non-nuU pruned Uee. Local areas for a two- 
30 dimensional input signal are shown in Figs. lOA and lOB to which reference is now briefly 
made. For signals upon which two-dimensional separable sampling was performed, the 
neighborhood is typically an NxN square of data sunounding the current node 35, where N 
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is typically 3 - 7. This is illustrated in Fig. lOA. If N is odd, then the node 35 sits in the 
veiy center of the neighboihood. Othenvise, it sits slightly to one side. If two^ensional 
quincunx subsampling perfonned, for even numbered levels, the local area is as shown 
in Fig. lOA. For odd numbered levels, the neighborhood is an NxN diamond, as shown in 
Fig. lOB. 

Referring back to Figs. 8 and 9, in step 82, the compaction unit 28 selects a 
first node 35 of the non-null pruned forest as the cunent "nucleus". In step 86 a zero- 
CTOSsing location for the cunent local area sunounding the cunent nucleus is computed. If 
the location of the zero-crossing happens to be at the cunent nucleus (step 88). then the 
cunent nucleus is "zero-crossed" and there are no more iterations. Otherwise, the cunent 
nucleus is set to be at the location of the zero-crossing and the process is repeated for the new 
local area around the cunent nucleus. The iterations continue until the cunent nucleus is al 
the location of a zero-crossing or a maximum number of iterations (typically 5 to 10) has 
been reached. The last value for the zero-crossing, no matter how far from the original node, 
is saved (step 94) and the process repeated for each node of the non-null pnmed forest 

The equations for detennining the location of the zero-crossing are provided 
heieinbelow in the section labeled "Zero Crossing Equations". 

The remaining steps of Fig. 8 create "zercncrossed" uees which contain only 
those nodes which are zero-crossed or which have a zero-crossing within their local area. 
Specifically, for each non-null pruned uee, the compaction unit 28 scans the first resolution 
signal HI. Whenever (step 102) it finds a "zero-crossed node" (i.e. one on which the zero- 
crossing occurs directly) fiom the list produced in step 90, it checks if the node is aheady a 
node in some zero-crossed tree (step 104) in which case, it stops the present detenninations 
and adds a temporary list STL to the zero-crossed tree to which the node belongs (step 1 10). 
If it is not already part of the zero-crossed tree, the unit 28 saves the node in the temporary 
list STL (step 106). 

Unit 28 then sets the cunent node to its parent node (step 108) and then has 
to determine if the parent node should be added to the zero-crossed tree. It repeats steps 104 
and 106 and, if the parent is not a root, the foDowing criteria are checked: (step 112) is the 
parent node zeio-crossed, (step 114) does it have a zero-crossing in its local area based on 
the infonnation from the method of Fig. 9; or (step 116) are its parent and its chfld both zero- 
crossed, b other words, step 116 indicates that the parent is a small gap in a branch of zero- 
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CTOSsed Dodes and is a lelaxation criterion vis-a-vis the other criteria. 

If the parent node fumils one of the criteria, it is added to the temporary list 
STL aBd the next parent is taken. The determination is repeated untfl a datapoint from the 
root signal is reached (step 1 18) or when a node which is already part of a zero-crossed tree 
is found. Otherwise, the determination for this node of the fust high frequency signal level 
HI is fmished and the next node is taken. When the nodes of HI are exhausted, the process 
is repeated for the other high frequency signals, starting from level H2. 

The nodes are scanned in quanta of local area size. Thus, the first resolution 
level HI is divided into units the size of one local area. For a two-dimensional input image, 
Oie units are squares or diamonds for the odd levels from quincunx subsampling. During a 
first pass, the daiapoinls in the same position in each local area, for example in the top left 
corner, are checked. In the next pass, the datapoints from the next position, for example, the 
next position in the top row. are checked. 

The daupoints are first checked that they are nodes of the non-null tree. If not, 

15 the next datapoint (i.e. from the next local area) is taken. 

The result of the operations of Fig. 8 are zero-crossed trees having therein only 
those datapoints which are at or near transients. However, a human typicaUy views a 
transient as having a certain thickness. Thus, the zero-crossed trees do not contain enough 
information. They must he filled with those nodes of the non-nuU pruned trees which are 
20 near the nodes of the zero-aossed trees. 

Therefore, in step 74 of Fig. 7, the compaction unit 28 adds data to the nodes 
of the zero-crossed trees. The operations performed in step 74 are listed, in flow-chart 
format, in Fig. 11. 

The nodes of the new trees, called herein "frilly pruned uees", are first set to 
25 be the nodes of the zero-crossed trees (step 130). 

The compaction unit 28 then scans the nodes N of each non-null pruned tree 
to detennine which nodes are either already part of the structure (step 132) or are part of the 
local area of at least one node of the corresponding zero-crossed tree (step 134). A list M 
is produced of the nodes of the zero-crossed uees whose local area the current non-null node 
30 N is in. 

If a node is neither part of the structure or part of the local area of a node, a 
new node is tested (step 136). If M has one member, m*, the cunent node N is added to the 
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fully pruned uee al node m* (step 138). If M has many members (le. the node N is part of 
many local aieas), the node m which is closest to the current node N is selected as m* (step 
140) and cuireut node N is added to the fully pruned tree at node m* (step 138). Hie process 
is then repeated for aD nodes in each non-null pruned tree. 

The process of Fig. 11 ensures that, if a node of a non-ntill tree should be 
included in the fully pruned tree, it is included at only one node. Thus, the fully pruned tree 
contains the areas of high information content in a non-redundant manner. 

The convener 40 (Fig. 3) of the decompression unit 22 reproduces the wavelet 
transform data W by placing the value of zero at every pixel not included in the list The 
inverse wavelet transformer 42 then reproduces the input signal from the uncorapacted 
wavelet transform data W . 



Nullilier Eouations; 
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W,(.l,ij>\WOm (3) 
W^(l,ij)=W,(l,ijym,ij) (4) 



Wmji^[5^.Q.5\ (5) 



where bcJ is the floor of X which is the maximum integer smaller than X. 
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(9) 



{WO,ij) otherwise 

Where mJJ) ^ a function of "iraponance" which may be fixed for each level 
( m,iJ>R(.^ ) or '^^ depeDdenl (for example, the human visual response funcUon may 
5 serve as the function R for image data). Q is a predetermined quantization step and W^^ 
is a predetermined maximum number of bins in the histogram. 



Zero Crossin p F-ouations: 

The foUowing is noted about the Kto crossing equations: a) they are writien for a two- 
dimensional input signal, b) as the zero crossing is done within the level, we drop the level 
index, c) W(ij) denotes a detailed signal daiapoint and d) the operations are done for a w 
wiihin a local area nucleus. 
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w>0 



w>0 



w>0 



w<0 



(14) 



w<0 
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w*0 



(15) 



(16) 



(17) 



ll will be appreciated by persons skilled in the art that the present invention 
is not limited to what has been particularly shown and described hereinabove. Rather the 
scope of the present invention is defined by the claims which foUow. 
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CLAIMS 



1. A system for compaction and reconstruction of wavelet transform data of an 

input signal, the wavelet transform data fonned of a plurality of resolution levels, the system 
5 comprising: 

a compaction unit, operative on forests of trees, each tree having branches 
which relate one datapoint in a low jesolution level with the entirety of its descendant 
datapoints in higher resolution levels, for pruning branches from said trees and trees from said 
forest to remove all but areas of high information content in said input signal and for 
10 producing thereby pruned trees; and 

a converter for convening said pruned uees to reconstructed wavelet transform 

data. 

2^ A system according to claim 1 wherein said compaction unit comprises a non- 

] 5 null pruning unit for removing branches of said trees having datapoints thereon whose wavelet 
transform values are below predetermined threshold levels and for producing thereby non-null 
pnined trees. 

3, A system according to claim 2 and wherein said compaction unit additionally 

20 comprises a zero-crossing pruning unit for determining locations of transients in said input 
signal and for creating zero-crossed pruned trees having therein only datapoints of said non- 
null pruned trees which are at or near said locations of transients. 

4 A system according to claim 3 wherein said datapoints of said zero-crossing 

25 pruned trees form nodes and wherein said compaction unit further comprises means for 
associating datapoints of said non-null pruned Uees with nodes of said zero-crossed pruned 
trees, wherein each datapoint can be within a local area of many nodes but is associated with 
only one of them. 

30 5, A system according to claim 2 and including a nullifier, operative prior to said 

compaction unit, for nulling those datapoints of said wavelet transform data which are below 
said predetermined threshold levels. 
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6. A system according to claim 3 and wherein said zero-aossed pnming unit 

includes means for determining locations of zeio^ossings in said wavelet transfonn data. 

5 7. A method of compaction and reconstruction of wavelet transform data of an 

input signal, the wavelet tiansform data formed of a plurality of resolution levels, the method 

comprising the steps of: 

pnming branches from Uees and trees from forests to remove aU but areas of 
high information content in said input signal and producing thereby pruned trees, wherein 
10 each tree has branches and each branch relates one daiapoint in a low resolution level with 
the entirety of its descendant datapoinis in higher resolution levels; and 

convening said pruned uees to reconstnicted wavelet Uansfoim data. 

8. A method according to claim 7 wherein said step of pruning comprises the step 
15 of removing those branches of said Uees having datapoints thereon whose wavelet Uansform 

values are below predetermined threshold levels and producing thereby non-null pruned uees 

9. A method according to claim 8 and also comprising the step of determining 
locations of uansients in said input signal and creating zero-crossed pruned trees having 

20 therein only datapoints of said non-null pruned uees which are at or near said locations of 
Uansients. 

10. A method according to claim 8 wherein said datapoints of said zero-crossed 
pruned trees form nodes and wherein said step of pruning further comprises the step of 

25 associating datapoints of said non-null pruned uees with nodes of said zero-crossed pruned 
uees. wherein each daiapoint can be within a local area of many nodes but is associated with 
only one of them. 

11. A method according to claim 8 and including the step, operative prior to said 
30 compaction unit, of nulling those datapoints of said wavelet Uansform data which are below 

said predetermined threshold levels. 
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12. A method according to claim 8 and wherein said step of determmmg 

the step of determining locations of zero-crossings in said wavelet transform data. 
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